APPENDIX B 


Effects of Aircraft Wake Dynamics on Measured and Simulated 
NO x and HO x Wake Chemistry 

D.C. Lewellen and W. S. Lewellen 


Journal of Geophysical Research 
Vol. 106, No. D21, pages 27,661-27,672 
16 November 2001 


12 



Effects of aircraft wake dynamics on measured and 
simulated NO x and HO x wake chemistry 

D. C. Lewellen and W. S. Lewellen 

Department of Mechanical and Aerospace Engineering, West Virginia University, Morgantown, 
West Virginia, USA 

Abstract. High-resolution numerical large-eddy simulations of the near 
wake of a B757 including simplified N0 X and HO x chemistry were performed 
to explore the effects of dynamics on chemistry in wakes of ages from a 
few seconds to several minutes. Dilution plays an important basic role in 
the NO x - O3 chemistry in the wake, while a more interesting interaction 
between the chemistry and dynamics occurs for the H0 X species. These 
simulation results are compared with published measurements of OH and 
HO2 within a B757 wake under cruise conditions in the upper troposphere 
taken during the Subsonic Aircraft Contrail and Cloud Effects Special Study 
(SUCCESS) mission in May 1996. The simulation provides a much finer 
grained representation of the chemistry and dynamics of the early wake 
than is possible from the 1 s data samples taken in situ. The comparison 
suggests that the previously reported discrepancy of up to a factor of 20 - 
50 between the SUCCESS measurements of the [H02]/[OH] ratio and that 
predicted by simplified theoretical computations is due to the combined 
effects of large mixing rates around the wake plume edges and averaging over 
volumes containing large species fluctuations. The results demonstrate the 
feasibility of using three-dimensional unsteady large-eddy simulations with 
coupled chemistry to study such phenomena. 

1. Introduction 

The possibility that future increases in the commer- 
cial aircraft fleet may have significant environmental 
and climatological effects is being widely studied [e.g., 

Fenner et a/., 1999; Schumann , 1998; Fnedl , 1997]. 

While the volume of aircraft emissions is not large com- 
pared to ground-based sources, these emissions are in- 
jected directly into the relatively stable upper tropo- 
sphere and lower stratosphere where small emissions 
may have potentially significant effects. A number of 
field programs have been conducted, e.g., Subsonic Air- 
craft Contrail and Cloud Effects Special Study (SUC- 
CESS) [Toon and Miake-Lye , 1998], SASS Ozone and 
Nitrogen Oxide Experiment (SONEX) [Singh et oi., 

1999], and AEROXOX [Schumann, 1997]. which use in 
situ measurements in wakes as a means to better deter- 
mine the aircraft emissions and understand the chemi- 
cal kinetics within the wake plume. Our subject in this 
work, the effects of wake dynamics on wake chemistry 


measurements, is an important source of uncertainty in 
interpreting the results of such field measurements or 
results from atmospheric chemistry models. When the 
sampling volumes being measured (or modeled) con- 
tain large fluctuations of some chemical species, the av- 
erage concentration measured and chemistry inferred 
may be quite different from the local results occurring 
within that volume. This is particularly true when re- 
action rates between relevant species are faster than 
turbulent mixing rates, allowing strong anticorrelated 
distributions to occur. This necessarily reduces the ef- 
fective reaction rate of average concentrations, making 
those rates dependent on mixing rates in the fluid [St/fces 
et ai., 1992; Schumann , 1989]. 

Here our focus is on the dynamics and chemistry 
within a wake plume from an age of a few seconds to 
several minutes. This is a regime which has been heavily 
studied in in situ measurement campaigns. It includes 
the partial capturing of the exhausts by the trailing 
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vortex pair behind the aircraft, the vortex-dominated 
dynamics which dramatically expand the wake plume 
both vertically and horizontally, and the subsequent 
buoyant oscillations of the plume after the vortex dy- 
namics subside. We do not consider the chemistry in 
the earliest jet exhaust phase where temperatures are 
initially very high but rapidly dilute toward ambient 
levels, or the ultimate fate of the wake plume and its 
chemical constituents at late times where the dynamics 
are dominated by ambient atmospheric turbulence, al- 
though both are likely to prove important in assessing 
the environmental effects of aviation. 

During this time period we are particularly inter- 
ested in the impact of species fluctuations on scales that 
are typically too small to be measured by trailing air- 
craft moving at speeds often in excess of 200 m s 
Accordingly, we utilize a three-dimensional (3-D), un- 
steady, high- resolution model of the aircraft wake which 
is capable of resolving dynamical and/or chemical fea- 
tures on scales of a few meters. The large-eddy simu- 
lation (LES) code we employ for the fluid dynamics is 
the same one used by Lewellen and Lewellen [1996] to 
present detailed dynamics of the vortex pair breakup in 
an aircraft wake; by Lewellen et al. [1998) to compare 
numerical simulation results of wake dispersion to lidar 
observations; and by Lewellen and Lewellen [2001] to 
examine effects of wake dynamics on contrail evolution. 
Given the computational requirements of the resolution 
we employ, we can include only a limited set of chemical 
reactions and species within our simulations. Nonethe- 
less, this proves sufficient to demonstrate some poten- 
tially important effects, while confirming the practical 
feasibility of coupling detailed dynamics and chemistry 
using LES. 

Of particular interest to us in this study were a set of 
HO* measurements taken during SUCCESS [Ton et al ., 
1998]. These showed an apparently anomalously high 
ratio of [H0 2 ]/[0H] within young wake plumes, 20 - 50 
times greater than expected by simple chemical compu- 
tations. The authors suggested that this effect might 
be due in part to the strong anticorrelation between 
NO and H0 2 in the turbulent fluctuations within the 
wake occurring during these measurements. To quan- 
titatively test this possibility, we have included HO x 
chemistry in addition to basic NO x - ozone chemistry 
in our reaction set for the present work, and have cho- 
sen to simulate a B757 wake with nominal conditions 
and emissions in agreement with those reported during 
the SUCCESS H0 X measurements. 

Our approach here is complementary to other recent 
studies of chemistry in aircraft wake plumes. These in- 


clude studies of detailed chemistry in the early hot jet 
regime coupled with simple dynamics [e.g., Mxoke-Lye 
et al, 1993; Menon and Wu , 1998; K archer et al. , 1998]; 
studies with detailed chemistry in the wake vortex 
regime with dynamics modeled in two dimensions [e.g., 
Anderson et al 1996]; studies of older plumes with de- 
tailed box model chemistry but only rudimentary dy- 
namics included via some prescribed dilution history 
[e.g., Karol et a/., 1997]; studies with detailed 3-D air- 
craft wake dynamics but only passive tracers for the en- 
gine exhausts [e.g., Lewellen and Lewellen , 1996; Gerz 
and Ehret , 1996; Gerz et al, 1998; Lewellen et al. , 1998]; 
and studies of ice phase “chemistry” (contrails) inter- 
acting with wake dynamics at many different levels of 
approximation [e.g., Gierens , 1996; Chlond , 1998; Suss - 
mann and Gierens , 1999; Gierens and Jensen , 1999; 
Lewellen and Lewellen, 2001; Sussmann and Gierens , 
2001]. Aircraft wake dynamics alone, ignoring the fate 
of the engine exhaust products, have also been a sig- 
nificant area for research because of potential effects on 
closely following aircraft and the resulting safety issues 
for takeoff and landing at airports. A recent review of 
these efforts, including some large eddy and direct nu- 
merical simulations, has been made by Spalart [1998]. 

In section 2 we provide an overview of our LES model 
and the conditions simulated. In section 3, we discuss 
the simulation results, beginning with the basic wake 
dynamics and NO x - ozone chemistry, as illustrated by 
the simulations, and ending with the dramatic effects of 
dynamics on measured HO x chemistry in the wake. A 
summary and some final comments follow in section 4. 

2, Model Overview and Simulation 
Conditions 

In order to simulate detailed wake dynamics, species 
dispersion, and the accompanying chemistry, we use 
the technique of numerical large-eddy simulation. The 
range of length scales involved in this problem is much 
too great to permit a direct numerical simulation of all 
scales with currently available computer resources. In- 
stead, we solve the Navier-Stokes equations directly for 
the most important scales in the How (from smaller than 
the wake vortex core size to many times the aircraft 
wingspan), but only model the effects of the smaller, un- 
resolved eddies. Our simulation is of a 13757 flying in the 
upper troposphere under cruise conditions, nominally 
that relevant for the SUCCESS observations. Details 
of our LES model applied to aircraft wake dynamics 
have been given by Lewellen and Lewellen [1996, 2001], 
along with a discussion of the many choices and com- 
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promises necessarily involved in efficiently tackling such 
simulations. A summary of the model and simulation 
conditions is given below, first for the fluid dynamics 
and then for the chemistry added specifically for the 
present work. 

2.1. Fluid Mechanics 

The LES code we employ for aircraft wake simula- 
tions is a modified version of one we use for bound- 
ary layer cloud modeling [e.g., Lewellen and Lewellen , 
1996], which in turn is based on the implementation of 
Sykes and Henn [1989]. It is a 3-D finite difference im- 
plementation of the incompressible Navier-Stokes equa- 
tions in the Boussinesq approximation, second-order 
accurate in space and time. It incorporates a piecewise 
parabolic model algorithm for the ad vert ion of thermal 
energy and species concentrations. The pressure field 
is obtained by directly solving the Poisson equation us- 
ing fast Fourier transforms in the downstream direction, 
the method of Famell [1980] for nonuniform grids in the 
cross-stream direction, and a tridiagonal solver in the 
vertical. 

The simulation of the interacting trailing vortex pair 
in an aircraft wake presents a particular challenge be- 
cause both small grid spacings and large domains are re- 
quired. A grid spacing of a few percent of the wingspan 
or less is required to resolve the trailing vortex cores 
and the transport of the engine exhausts around and 
into them. At the same time, line vortices have a sizable 
long-range influence on the flow field, so domain widths 
of several wingspans or more are required to minimize 
unwanted boundary interactions. We employ three de- 
vices to help to meet these conflicting demands: (1) a 
stretched grid spacing in the vertical and cross-stream 
directions to achieve fine spacing where it is needed 
within the wake while maintaining large domain sizes; 
(2) vertical tracking of the wake to insure that the fine 
grid portion of the domain encompasses the descending 
wake vortices; and (3) different grids (six in the case 
presented here) for different segments of the simulations 
as the wake evolves and the optimum grid requirements 
change. A variable time step was used in the simula- 
tion, with values ranging from 0.003 to 0.02 to 1 s for, 
respectively, the early jet, early vortex, and late time 
stages discussed below. 

We employ a subgrid model which uses a quasi- 
equilibrium, second-order turbulence closure scheme 
with the subgrid turbulent kinetic energy carried as a 
dynamical variable and the maximum subgrid turbu- 
lence length scale related to the numerical grid length. 
The damping of subgrid turbulence by a stable tem- 
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perature gradient or a stable rotational velocity field 
is approximately included via appropriate reductions of 
the subgrid turbulence length scale. The incorporation 
of rotational damping effects in the subgrid model is 
somewhat unusual, but useful for simulations of strong 
vortices in order to avoid unphysically large rates of core 
diffusion. The method is described in detail by Lewellen 
et ai . [2000], who used it extensively in tornado simula- 
tions. It was not included in our earliest aircraft wake 
simulations [Lewellen and Lewellen , 1996], but has been 
used since [Lewellen et al , 1998; Lewellen and Lewellen , 
20011. 

We employ periodic boundary conditions in all three 
directions. Effectively, this means there are neighboring 
wakes located at twice the distance to the boundary in 
the vertical and crosswise directions evolving in parallel 
with our simulated wake. The stretched grid allows us 
to place the boundaries at sufficient distance to make 
the effects of these image vortices negligible. Period- 
icity in the vertical is modified by a temperature shift 
between the top and bottom of the domain to permit a 
mean stratification to be imposed. In imposing periodic 
conditions in the downstream direction we are treating 
the wake portion within our domain as if it were laid 
down instantaneously. This is an appropriate approxi- 
mation here because the aircraft flight speed is much 
greater than the velocity scales of the wake dynam- 
ics over the period treated (which excludes the region 
within a few wingspans of the aircraft.). The approxi- 
mate periodicity of the wake downstream is amply con- 
firmed by observations of contrail evolution from ages 
of a few seconds to several minutes. 

The initial conditions were chosen consistent with 
cruise conditions for a B757 flying at 11 km altitude at 
210 ms" 1 , with a mass of 87,000 kg (which is midway 
between the minimum and maximum masses for this 
aircraft). For an effective vortex pair separation of 29.4 
m this implies a circulation about the vortices of 380 
m 2 s 1 . The heat produced by the engines per meter 
of flight path was taken as 255 K m 2 , computed from a 
total fuel flow rate of 3.1 kg kin' 1 (in agreement with 
that quoted for the SUCCESS measurements by Tan 
et ai [1998]), a specific combustion heat of 13 M.I kg 1 
appropriate for kerosene- type jet fuels, and an assumed 
jet engine efficiency of 30%. The integrated momentum 
in the exhaust jets per meter of flight path was chosen 
as 502 m 3 s l , computed from the assumed fuel flow 
rate, flight speed, and a specific fuel consumption of 
16.9 mg N -1 s'* 1 . The detailed cross-stream distribu- 
tions of the velocity fields and exhaust products depend 
on the roll up of the vorticity shed from the wing and 
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were taken from a UNTVVAKE simulation [Quackenbush 
et cd., 1996], with values scaled so as to match the de- 
sired circulation, jet momentum, engine heat, etc. 

As we will discuss in detail below, the fluctuations in 
exhaust distribution play an important role in affecting 
measured chemistry in the wake. Accordingly, we have 
at least approximately included the effects of exhaust 
jet turbulence into the initial conditions used for our 
chemistry simulation. This was done by performing an 
initial simulation of the early w’ake (without chemistry) 
from an age of about 0.5 s to 3.5 s, with high spatial 
resolution (0.33 m finest grid spacing), for a short (42 
m) downstream portion of wake, to develop some jet 
turbulence and smaller-scale fluctuations in the exhaust 
distribution. These results were periodically continued 
to provide the initial conditions for the full downstream 
domain employed for the chemistry simulation. 

To complete the initial conditions for our simu- 
lations, the aircraft wake is superimposed with an 
ambient atmospheric field (the result of a separate 
large-eddv simulation), as described by Lewellen and 
Lewellen (1996]. The field was generated in the desired 
domain, with a nominal value of potential temperature 
stratification for the upper troposphere (2.5 K km -1 ) 
and no mean wind shear. Large initial velocity and 
temperature perturbations were applied and the subse- 
quent evolution followed via LES until the turbulence 
and gravity waves had a chance to adjust to the stratifi- 
cation, and their amplitude had decayed into a reason- 
able range (turbulent dissipation rate of order 2 x 10“ 5 
m 2 s“ 3 ). 

Table 1 summarizes the grids employed for the pri- 
mary simulation. The choice of grid for the given seg- 
ment of the simulation was correlated with the wake 
vortex dynamics then occurring, in order to best appor- 
tion the limited number of grid points available given 
the physics that had to be resolved. As the dynam- 
ics progress the wake expands and turbulence levels 
fall, so that the finest grid spacing can be increased on 
subsequent grids but must cover larger portions of the 
domain. Accordingly, when the wake dynamics would 
threaten to outgrow the fine-grid region of the domain, 
we simply interpolated the species and flow fields to a 
more suitable grid. The domain size was doubled for 
the fourth grid. To do this, the final field generated 
on the third grid w’as periodically doubled in the down- 
stream direction and added to an ambient turbulence 
field generated by a separate LES on the doubled do- 
main. 

As noted previously, the vertical and cross-stream 
domain edges were chosen far from the wake to mini- 


mize boundary effects and allow a more realistic ambi- 
ent turbulence field. There is a different constraint on 
the downstream domain size as discussed at length by 
Lewellen and Lewellen [1996]. The downstream domain 
size mast be commensurate with the period of the basic 
wake decay mode, the Crow instability (discussed be- 
low). For a given aircraft, there is considerable variabil- 
ity in the observed periodicity, depending on the wave 
and turbulence spectrum in the ambient atmosphere. It 
is not unusual to see a factor of 2 variation in period 
within a single contrail in the sky. Nonetheless, there 
is a preferred range of wavelengths encountered. Here 
we choose a downstream domain dimension of 210 m, 
expecting a Crow period of approximately that length, 
which for the B757 is a little less than the wavelength of 
maximum amplification deduced by Crow [1970] from 
a linear stability analysis, and a little greater than the 
most favored range found in the direct numerical sim- 
ulation study of Spalart and Wray [1996]. In previous 
work we have found little variation of the basic wake 
decay (rate of wake dispersion, overall vertical and hor- 
izontal extent, wake lifetime, etc.) with downstream 
domain size. The interaction between a given Crow pe- 
riod and its neighbors is only weakly affected by the 
detailed structure of those neighbors, so replacing them 
with periodic images proves to be a reasonable approx- 
imation. Providing for multiple Crow periods to be 
simulated does give additional turbulent realizations of 
these structures. We performed an additional simula- 
tion analogous to the one presented here but with a dif- 
ferent ambient turbulence field (with lower turbulence 
intensity) to estimate the variation in our chemistry 
results between different turbulent realizations of the 
same basic process. For the summary results presented 
below, the two sets of results were all very close. 

2.2. Chemistry 

Given the number of grid points in our simulations, 
we use a relatively simple (and therefore numerically 
less costly) chemical reaction set, concentrating on im- 
portant reactions within the wake plume for ages from 
a few seconds to a few minutes. The chemical species 
considered are NO. NO2, HONO, HNO3, O3, OH, HO2, 
CO, and CH*, with basic reactions: 


NO + 0 3 

— 4 N0 2 + 0-2, 

(1) 

NO2 4 - 4- ()■> 

— > NO + O,, 

(2) 

HONO + h u 

— ► OH 4- NO, 

(3) 

OH + NO 

— ► HONO, 

(4) 

OH + NO* 

—4 HNOj, 

(0) 
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Table 1. Grids Used for the B757 Chemistry Simulation 0 


Grid 

Domain, m 

Number of Grid Points 

Fine Spacing, m 

Time period, 

0 

660x660x42 

162x154x130 

0.33x0.33x0.33 

0.5-3. 5 

1 

660x800x210 

124x116x130 

0.8x0.8xl.6 

3.5-50. 

2 

660x800x210 

120x121x162 

l.xl.xl.3 

50-100 

3 

660x900x210 

152x134x130 

1.6xl.6xl.6 

100-240 

4 

1200x1300x420 

100x139x122 

3.3x3.3x3.5 

240-500 

5 

1200x1300x420 

94x119x66 

6.5x6.5x6.5 

500-900 


“Domain dimensions, number of grid points, and fine grid spacings are given for the (cross- 
stream)x (vertical)x(downstream) directions, respectively. Grid 0 was used for an initial simula- 
tion without chemistry to generate initial conditions including turbulent jet fluctuations for the 
main simulation. 


H 0 2 + NO 

— > OH + NO2, 

(6) 

OH + CO + 0 2 

— * CO2 + IIO2, 

( 7 ) 

OH + 0 , 

— ► HO2 + 0 2 , 

(8) 

H 0 2 + Os 

— -»• OH + 2O2, 

( 9 ) 

OH + CH4 +... 

— ► HO2 + NO2 + ... 

(10) 


The reaction rates employed (&no+o 3 = 2.39x10 
fcoH+NO = 4.63x10 &oh+no 2 = 9.18x10 2 , 

^HO^+NO — 8.4xl0~ 2 , fcoH+CO = 1.29x10 ^OH+Oa 
= 1.59xl0“ 4 , fcHo 2 -f0 3 = 8.32xl0~ fi , fcon+CH* — 
5.2xl0 -6 , all in (ppb s) _l ) were taken from DeMore 
et al. [1997] for conditions of 217 K and 0.224 atm, ex- 
cept for the OH 4- NO*> reaction rate, which was taken 
from Brown et al. [1999]. Nominal daytime values at 11 
km altitude were chosen for the photolysis rates ( Jno? 
= 0.01 s' 1 and Jhono = 0.0033 s" 1 ). Reaction (10) is 
an effective one for the cycling of OH to HO2 through 
the intermediary CH3O2 in the presence of NO. It is 
only of minor importance to the HO x chemistry in the 
regime being st udied (being dominated by reaction (7)), 
but is included approximately at little extra numerical 
cost by taking the CH 4 concentration as constant. The 
variations of the rate constants due to temperature fluc- 
tuations within the wake plume have not been included 
here. For the wake ages considered, the variations in 
temperature encountered are generally of order a de- 
gree or less, giving changes in rate constants of a few 
percent or less (which is generally below the level of 
uncertainty in the rate formulas themselves). 

In choosing the initial chemical concentrations, we 
did not try to match a particular SUCCESS flight con- 
dition, but did pick values within the ranges measured 


during the SUCCESS flights: ambient background lev- 
els of HO2 between 10-20 pptv and OH between 0.3-0. 5 
pptv [Tan et al y 1998]; ambient levels of NO* = 70 
pptv, CO = 120 ppbv, 0 3 = 65 ppbv, and CH 4 = 1750 
ppbv [Jaegle et al 1998]- The actual values used for 
the HONO, OH, and H0 2 concentrations (HONO = 
0.36 pptv, OH = 0.42 pptv, and H0 2 = 15 pptv) and 
the split of NO x into NO and NO2 (NO = 55 pptv and 
NO2 = 15 pptv), were also chosen so that the back- 
ground concentrations of any species did not drift sig- 
nificantly over the course of the simulation (realizing 
that other reactions besides those explicitly treated are 
important in setting the background equilibrium con- 
centration over timescales much longer than we consider 
here). 

The initial concentration used for each species con- 
sisted of a background level plus a plume concentration, 
with the latter chosen consistent with measured emis- 
sion indices, fuel consumption rates, and flight speeds 
for the SUCCESS flights. In particular, the NO emis- 
sion index is taken as 7.5 g N0 2 kg -1 fuel [Campos 
et aL, 1998], with a total fuel usage per flight path of 
3.1 g m -1 . The volume ratio of CO above background 
in the exhaust to NO x was set at 0.306, consistent with 
the median emission index (1.55 g kg J ) measured on 
May 4 during SUCCESS [ Vay et af., 1998]. The ratios 
of HONO and N0 2 to NO x in the jet are set at nominal 
values of 0.04 and 0.1, respectively. The initial plume 
excesses over background for HNO3, O3, CH 4 , OH, and 
HO 2 were taken to be zero. 

The biggest numerical costs of including chemical re- 
auctions within our LES are the added memory required 
for the 3-D species fields and the added computation 



6 


LEWELLEN AND LEWELLEN 


involved in adverting and diffusing them as governed 
by the flow field. Accordingly, it is useful to minimize 
the number of species which must be carried as dy- 
namical variables. We reduce this number by two by 
identifying three combinations which are conserved by 
equations (1)— (10) (so that their plume concentrations 
above background always remain proportional to each 
other): NO + NO 2 - OH - HO 2 (essentially NO x within 
the plume); H0 2 + OH + HONO + HNO 3 (which we 
use after the computation to solve for HNO 3 ); and CO 
(which is conserved to a high degree of accuracy because 
the CO levels everywhere are so much greater than the 
OH levels which erode it via (7)). 

In numerically implementing (1)-(10), one point re- 
quires particular care. Some of the reactions involv- 
ing OH and H0 2 (equations (4)-(6)) proceed on a very 
short timescale within the core of the wake plume where 
the NOx concentrations are high, shorter than the time 
step required to accurately follow the flow dynamics. 
We cannot simply use equilibrium values for OH and 
HO 2 because the timescale for the relevant reactions 
does not remain small (compared to the time step ap- 
propriate for the flow dynamics) within the edges of 
the plume. Instead, we update the OH and HO 2 chem- 
istry semi-implicitly. In short domain tests we have con- 
firmed that the simulated OH and 1102 concentrations 
evolve correctly in time within the plume edges, while 
giving the correct equilibrium values within the plume 
core, regardless of the time step chosen. The remaining 
reactions progress on timescales which are comparable 
to or slower than the flow timescales; consequently, the 
chemistry for the species other than OH and HO 2 are 
updated explicitly in time. 


3. Results 

The chief impact of wake dynamics on the chemistry 
is through plume mixing. This reduces concentrations 
of exhaust constituents, mixes fresh ambient air with 
plume air. and produces fluctuations in plume distribu- 
tions. It is appropriate, then, in explaining the simula- 
tion results, that we begin with a summary of a typi- 
cal wake decay history and its effects on passive plume 
dispersal. Simple NO* 0 3 chemistry is considered 
next before turning to HO x chemistry, where we find 
the most interesting effects of dynamics on chemistry 
illustrated by our simulations. 


3.1. Review of Basic Wake Decay Mode and 
Species Dispersion 

Figure 1 illustrates some of the simulated wake dy- 
namics out to an age of 10 min, as visualized by the 
cross-stream integrated NO x distribution in the wake 
plume. Since NO x is essentially a conserved species, the 
distribution shown is a result of dispersion only. The 
format is as described by Lewellen et al . [1998], inspired 
by scanning lidar measurements of wakes. The wake is 
sampled as if it were being advected at a steady rate 
bv a mean wind aligned with the wake. The horizontal 
axis then varies over time and downstream distance so 
that both the temporal evolution and spatial structure 
can be seen. The “drift velocity” determining the unit 
of distance per unit of time has been chosen to be 9.5 m 
s~ 1 , which is in a convenient range for displaying both 
the temporal and spatial structure. In the figure the 
vertical distance scale has been expanded by a factor of 
3 relative to the downstream distance scale for clarity. 

The main features of the wake evolution can be fol- 
lowed in this figure. We give only a brief summary of 
the dynamics here; more detail may be found in our pre- 
vious papers. The early wake dynamics are dominated 
by the vortex pair which results from the lift gener- 
ated by the aircraft's wings. For the first minute or so 
the vortex pair descends rapidly. Initially, the engine 
exhaust distribution is stretched and advected by the 
velocity field accompanying the vortex pair, giving a pe- 
riod of very rapid dilution. Subsequently, the exhausts 
are partitioned into three different populations: rela- 
tively high concentrations in and immediately around 
the vortex cores where mixing rates are low; midlevel 
concentrations and mixing rates within the ellipse of 
fluid surrounding and descending with the vortex pair; 
and lower concentrations and higher mixing rates within 
a fraction of the exhausts which is not ultimately en- 
trained into the falling vortex system (or becomes de- 
trained from it) and remains close to flight level. The 
fraction of exhausts finding its way into each population 
is most strongly influenced by the engine placement of 
the particular aircraft, but is also influenced to a lesser 
extent by the atmospheric turbulence levels, shear, and 
stratification. 

While the vortices descend, perturbations of the vor- 
tices from the ambient atmosphere (and/or unsteady 
wing loading or jet turbulence) grow via the Crow in- 
stability [Crow, 1970] . The rate of growth depends 
primarily on the wavelength of the perturbation, the 
vortex spacing, and circulation. In a turbulent atmo- 
sphere the vortices see a superposition of perturbations, 
and the instabilities for different wavelengths compete 
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Figure 1. Downstream space/time versus height plot of cross-stream integrated plume NO* (in ppm m) for a 
simulated B757 wake. 


with each other. Those which are in the range of most 
rapid growth and have the largest initial perturbations 
eventually win out. 

The growth of the Crow instability proceeds until the 
vortices touch, reconnect, and form elongated vortex 
rings. This occurs at a wake age of around 80 s in the 
present B757 simulation. During the reconnection pro- 
cess some of the exhausts escape from the vortex system 
and buoyantly rise. The distribution of the remainder 
is significantly scrambled so that the peak concentra- 
tions no longer provide a clear signature of the vortex 
dynamics. The vortex rings continue to descend and, at 
the same time, open up in the cross-stream direction, 
thereby significantly increasing both the vertical and 
horizontal extents of the wake. They can further inter- 
act with themselves, even sometimes relinking again to 
form two rings from one. This, together with the contin- 
ual diffusion of the vortex cores (due to the smaller- scale 
atmospheric turbulence and interaction with the ambi- 
ent stratification), eventually leads to the demise of the 
organized vortex system. In the present simulation the 
descent of the vortex rings (or their remnants) continues 
until an age of approximately 230 s, with about half of 
the total wake descent (of ~320 m) occurring after the 
pair initially relinks. It is the descending vortex rings 
which (for appropriate ambient humidity) give rise to 
the periodic series of puffs often seen in contrails. 

The coherent w T ake dynamics do not end with the 
dissolution of the wake vortices. The positive buoy- 
ancy acquired by the vortex system’s descent through 
the stratified atmosphere (together with the heat in the 
initial engine exhaust) drives a strong Brunt-Vaisala 
oscillation. In the absence of significant wind shear, 


we find that this oscillation drives the dominant turbu- 
lence that mixes the wake plume until it damps, after 
1-3 Brunt-Vaisala periods, to the level of the ambient 
atmospheric turbulence. The Brunt-Vaisala oscillation 
can be seen in the long- timescale (10 min period) ver- 
tical modulation evident in Figure 1. 

There are variations in the details of this basic wake 
decay mode depending on aircraft and atmospheric con- 
ditions, but that illustrated is typical. We have found in 
simulations of wakes and observations of contrails that 
the dynamics of the Crow instability primarily govern 
the wake decay for a wide range of conditions in the 
absence of large wind shear. 

3.2. NO* - Ozone Chemistry 

Figures 2-4 summarize the evolution of several chem- 
ical species during the wake decay. All of the concen- 
trations plotted are excesses above (or below) the con- 
stant background level. Consider first the evolution of 
the NO, NO;, and IIONO concentrations. These are 
dominated by plume dispersion. The local peak con- 
centrations of these species fall sharply (Figure 2) for 
the first -25 s as the jet exhausts mix throughout the el- 
lipse of fluid falling with the vortex pair. Subsequently, 
they fall nearly linearly while the vortex-driven turbu- 
lence dominates and more slowly when the vortex sys- 
tem dissolves and the weaker Brunt-Vaisala sloshings 
take over. There is a noticeable enhancement in the di- 
lution rate of the peak concentrations for a period after 
vortex linking occurs, with its accompanying scrambling 
of the exhausts. The peak local concentrations are lo- 
cated within the lowest (vortex-dominated) part of the 
wake plume until about the 250 s point, after which the 
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Figure 2. Peak concentrations of chemical species in 
the wake relative to background versus time. 



Figure 3. Peak downstream averaged wake plume con- 
centrations relative to background versus time. 

peak is located near the original flight level, within the 
upper tvake (which, being less turbulent, dilutes at a 
slower rate). 

In an in situ measurement, a following aircraft ob- 
tains not a strictly local concentration, but the average 
over some finite sampling volume. A corresponding plot, 
of peak measured concentration versus plume age would 



Figure 4, Integrated wake plume concentrations above 
background per meter of flight path versus time. 

look more like that in Figure 3 (corresponding to 1 s 
samples of a chase aircraft maintaining a constant sep- 
aration behind the lead aircraft). The NO, NO2, and 
HONO concentration histories in Figure 3 are well cor- 
related with those in Figure 2, but with some noticeable 
differences. The downstream average peak concentra- 
tions drop much more steeply between about 50 and 
100 s as the vortices are twisted by the growing Crow 
instability. The peak measured concentrations remain 
located with the vortices during this period, but sub- 
sequently shift to the exhausts released by the vortex 
system during the linking process. These buoyantly rise 
to rejoin the upper wake at about 180 s. 

The plume-integrated amounts of NO, N0 2 , and 
IIONO vary much more slowly, in response to the wake 
chemistry itself (Figure 4). The plume NO is slowly con- 
verted to N0 2 via (1); the initial ratio of [N0 2 ]/[N0] 
= 0.11 will eventually rise to approach its background 
equilibrium value for the present conditions of 0.27 as 
the plume dilutes. The total HONO is slowly depleted 
via photolysis, with a corresponding nearly linear rise 
in HNO :i to approximately 2.6% of the plume NO* at 
000 s. The time evolution of the peak HNO3 concentra- 
tion (Figure 2) represents a competition between this 
production and the plume dilution. 

The OH and H0 2 behavior in the plume are quali- 
tatively different. As we will discuss in more detail be- 
low, the differences from background for these species 
axe more nearly constant within the plume, so that the 
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integrated contribution scales approximately with the 
plume volume. There is a sharp initial increase associ- 
ated with the early mixing of the exhaust jets by the 
vortices, a nearly linear increase for the next 200 s as the 
descending vortices and oscillating vortex rings stretch 
the plume volume both vertically and horizontally, and 
a slower rate of increase as the vortices dissolve and the 
Brunt- Vaisala oscillations take over and slowly subside 
(Figure 4) . 

The O3 evolution is more interesting, and we might 
expect more subtle interplay between the chemistry and 
dynamics for this species since both the NO photolysis 
and NO + O3 reaction (for early plume NO concentra- 
tions) occur with characteristic timescales of order 100 
s, a timescale which is also of importance for the wake 
dynamics as seen in Figure 1 . The destruction of O3 by 
NO within the plume is larger than the production via 
NO photolysis for the first ~170 s (Figure 4) and re- 
versed afterward so that an O3 deficit forms but is later 
replaced by an 0 3 surplus, with the age of the crossover 
point strongly influenced by the plume dynamics. The 
magnitudes above or below background are all small, 
however. At the given NO concentrations, equation (1) 
has difficulty competing with the level of plume dilution 
continually mixing in fresh O3 from the air outside the 
plume. If either the NO concentration had been an or- 
der of magnitude larger, or the plume dilution rates that 
much smaller, the chemistry and dilution rates would 
have been comparable and the interaction between the 
two would have been more complex. An even larger 
change would make the chemistry fast compared to the 
mixing, and a much stronger anticorrelation between 
NO and O3 would have developed(as seen in the works 
by Schumann [19891 and Sykes et al. [1992]). 

3.3. HO* Chemistry 

The effects of the three-dimensional wake dynamics 
are much more dramatic for the OH and HO2 chemistry, 
and can explain a puzzle found in the in situ observa- 
tions of Tan et al. [1998]. In a series of observations 
during SUCCESS. OH. H0 2 , and NO concentrations 
were measured within the wake plume of a B757 for 
wake ages from about 10 to 100 s. The measuring DC- 
8 also sampled its own wake at much later ages (about 
1200 s). As expected (and as we find here), the mea- 
sured OH levels sharply rose upon entering the plume, 
while the H0 2 levels sharply fell (rf. their Figure 1). 

The measured ratio of HO2/OH concentrations as a 
function of NO, however, was about 20-50 times larger 
than expected (cf. their Figure 2). As the authors 
argued, within the plume where the NO levels are high 
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the sources and sinks of OH and H0 2 are slow compared 
to the exchange reactions betw r een them, so that the 
ratio is, to good approximation, 

[H0 2 ] _ fcpH+co [CO] 4- frpH+cm [CH4 ] A fepH+Os [O3] 
[OH] fcHOa + No[NO] 4- fcHOa+OaiOs] 

( 11 ) 

It is this prediction which is sometimes greatly exceeded 
by the observations. 



1 NO ) ( ppbv) 

Figure 5. Apparent [H0 2 ]/[0II] ratio deduced from 
simulation data using 210 m sample lengths, versus the 
equilibrium expectation (thick line). Small points are 
sampled from wake plumes of ages 6 and 10 s; open 
circles for plume at 50 s. 

Figures 5 and 6 show our simulated [H0 2 ]/[OH] ratio 
for sample measurements. The format is as in Figure 2 
of Tan et al. [1998] (except for the log scale and sorting 
the data by wake plume age), with the data sampled 
in analogous fashion. The time resolution for the mea- 
surements was 1 s, which at the average flight speed 
of 210 m s 1 of the measuring aircraft means an av- 
erage over 210 m for each data sample. Accordingly, 
the simulated measurements in Figures 5 and 6 are 210 
m averages sampled along lines parallel with the flight 
track at many positions within the simulated wake, at 
ages as noted. As in the observations, the analogous 
[IIO2] / [OTT] measurement from the simulation shows 
values ranging up to 50 times greater than the local 
equilibrium predictions (thick line). 

In our simulated chemistry there are two main effects 
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Figure 6. As in Figure 5, but for wake plume ages of 
100 s (dots), and 250 s (circles). 


responsible for this discrepancy, both of which critically 
involve the wake dynamics. The first stems from the 
difficulty with using spatially averaged samples in (11) 
when there are significant species fluctuations within 
the sampling volume. The second arises because the 
equilibrium assumption in (11) is no longer valid around 
the edges of the plume. Let us consider each of these 
effects, in turn, in more detail. 

Figure 7 shows the evolution of the simulated OH 
and HO 2 concentrations on a vertical strearnwise slice 
through the wake. The format is as in Figure 1 ex- 
cept the fields are not cross-stream integrated but lie 
on a single plane parallel to the wake, initially situated 
midway between the aircraft centerline and one of the 
trailing vortices. The vertical and downstream lengths 
are plotted at the same scale in this figure. In most of 
the plume volume at any given time OH is nearly con- 
stant. II 0 2 is very small, and (though it is not readily 
apparent from the figure) equation (11) is well satisfied. 
However, it is clear that our 210 m downstream samples 
will in general include a wide range of plume conditions, 
sometimes both high plume concentrations and air out- 
side of the plume. This sampling effect influences the 
measured HO 2 in the plume more than the measured 
OH, and the ratio still more because OH and H0 2 are 
strongly anticorrelated. This can be seen as follows. 
Within the plume, where NO is large and the OH H0 2 


exchange reactions are fast, the dominant source of OH 
and HO 2 is HONO photolysis, and their concentrations 
are well represented by the equilibrium values: 


[OH] 

[H0 2 ] 


«/hono [HONO] 


fcoH+No[NC)] 4* fe3H+N0 2 [N0 2 ] 

fOTTl fo H +C Q t C °] 

1 J fcHo,+No[NO ] 5 


, ( 12 ) 
(13) 


where the last equation follows from keeping the domi- 
nant terms in (11). The HONO, NO, and NO 2 concen- 
trations vary strongly within the plume, but this is due 
mostly to dilution, which acts on each of these species 
in essentially the same way. The ratios of these species 
(along with OH via (12)j then vary little within the 
plume and relatively slowly in time. In contrast, HO 2 
in the plume varies approximately as [OH]/[NO] since 
CO is approximately constant (there is a plume contri- 
bution to CO but it is generally well below the back- 
ground contribution except for the peak values at the 
earliest times in our simulation, when they are com- 
parable). In the plume, then, we have approximately 
[H0 2 ]/[OH] - l/[NO] (which is why the prediction from 
(11) gives a nearly linear curve on the log-log plots of 
Figures 0 and 6). The averaging from the finite sam- 
ple size affects the measurements on the two axes in 
Figures 5 and 6 very differently. The correlation of the 
averages (l/[NO]> and ([NO]) differs greatly from the 
correlation of point samples of l/[NO] and [NO] be- 
cause of the nonlinearity of the function being averaged. 
When averaging over sample volumes containing large 
fluctuations of the NO field, l/([NO]) and (l/[NO]) are 
clearly significantly different functions, the latter be- 
ing disproportionally enhanced relative to the former 
by the contribution of regions with low NO. Moreover, 
if some air outside of the plume (where (12) and (13) 
do not hold) is included within the sampling volume, 
then <[H0 2 ])/([OH]) will be further enhanced because 
[H0 2 ]/[OH] in the background is so large (~35). It is 
for these reasons that this sampling effect leads to an 
increase in the measured ratio. 

The effects of averaging strongly fluctuating concen- 
trations over finite volumes will be reduced as the sizes 
of the sampling volumes are reduced, and removed en- 
tirely for point samples. Figures 8 and 9 are the same as 
Figures 5 and 6 but with the data now sampled strictly 
locally. This measurement is in much better agreement 
with the equilibrium predictions, showing how much the 
spatial averaging skewed the previous results. There 
are still significant differences, however, due to the in- 
adequacy of the equilibrium approximations in large 
regions around the plume edges (as can be seen from 
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Figure 7. Concentration fields relative to background on a downstream space/time versus height slice for the 
simulated B757 wake: (a) OH (in ppt) and (b) HO* (in -5 ppt units). 


the regions of enhanced OH in Figure 7). The mixing 
into the plume of background HO -2 can be an important 
source, via reaction (6), of OH in the plume (in addition 
to HONO photolysis as assumed in (12)). Here as well, 
where NO is lower, the exchange reactions between OH 
and HO 2 are not necessarily fast compared to the dy- 
namics. Again this effect leads to larger [H02]/[0H] 
ratio at a given NO level than predicted in equilibrium, 
because reaction (G) converting HO 2 to OH (pushing 
them to their equilibrium levels) is slowed relative to 
the dynamics mixing HO 2 into the region. 

The nonequilibrium effect is largest for the young 
plume where the vortex-driven mixing is most vigorous, 
and becomes steadily less important in the aging plume 
as the turbulence levels fall. In contrast, the effect from 
sampling finite volumes involves fluctuations on differ- 
ent scales at different times during the evolution, and 
remains important out to later times. At early times, 
which tend to cover larger values of NO. the sharp edges 
of the engine exhaust and its instabilities contribute 
most to the effect, with fluctuations on scales of order 


10 m or less. Subsequently, the engine exhausts get- 
smoothed out somewhat by diffusion as they are more 
completely roiled into the vortex pair. The dominant 
downstream species fluctuations here are of intermedi- 
ate scale arising from small axial pressure variations 
within the vortex cores (due to ambient turbulence). 
These fluctuations can be large for larger ambient tur- 
bulence levels (sometimes misrepresented as being due 
to vortex “breakdown” or “bursting”) but are modest 
here. We see reduced sampling effects in Figure 5 for 
the wake age of 50 s. From about 60 s and later, the dis- 
tribution is most strongly affected by the growth of the 
Crow instability of the vortex pair. The sampling effect 
in Figure 6 from these fluctuations (and presumably in 
the SUCCESS observations) is large because the sample 
length is approximately equal to the most favored wave- 
length for the Crow instability for this aircraft. At later 
times these induced fluctuations are slowly eroded down 
to a level consistent with the ambient turbulence as the 
wake slowly dilutes. By reducing the sampling size well 
below the Crow wavelength, the sampling effects could 
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Figure 8. As in Figure 5, but using local point samples. 



Figure 9. As in Figure 6, but using local point samples. 

be significantly reduced for all but the earliest wake 
times. 

It should be noted that we are almost certainly un- 
derestimating the effects of sampling over fluctuations 
in the very young wake versus those encountered in 


the SUCCESS measurements, for at least two reasons. 
First, we do not begin the simulation of the exhaust 
jets at their earliest and most intense phase, and so 
we underpredict the fluctuations from jet turbulence. 
Second, and probably more importantly, we have per- 
fectly aligned our sampling lines along the initial wake 
axis, while the DC-8 sampling behind the B757 during 
SUCCESS would have had difficulty remaining centered 
within the small cross section of the young wake plume, 
especially given the large wake turbulence present there. 

Finally, one apparent difference between Figure 2 of 
Tan et al [1998] and our results deserves comment. As 
we have noted, when our measured [H 02 ]/[ 0 H] ratio 
as a function of NO differs appreciably from the pre- 
diction of (11), it exceeds it (as is appropriate for the 
explanations of this effect given above). This is true for 
the bulk of the data points given by Tan et al. [1998] 
as well, but they also find a small fraction scattered 
below the expected values. These may just be stray 
data points (they quote a i-a precision for their data 
of 0.23 pptv for OH and 0.75 pptv for HO 2 ), indicate 
fluctuations not considered here in background levels of 
some species, or be due to the fact that the NO and 
HO* probes on the DC-8 were not colocated but on the 
starboard and bottom of the fuselage, respectively. 

4. Concluding Remarks 

This LES of the near wake of a B 75 7 provides a fine- 
grained chemical-dynamical representation of simplified 
NO x -HOx chemistry in the early wake. It demonstrates 
the feasibility of such simulations to investigate the po- 
tential impact of turbulent dynamic interactions with 
the wake chemical kinetics. By sampling the simulated 
data in a manner similar to that of in situ aircraft mea- 
surements, it is possible to interpret results obtained at 
the typical sampling rate of 1 Hz. In particular, com- 
parison of these simulation results with actual in situ 
flight measurements of OH and IIO 2 performed under 
cruise conditions in the upper troposphere during SUC- 
CESS in May 1996 provides a likely explanation for the 
apparent large discrepancy between measurement and 
theoretical expectation found by Tan et al. [1998]. The 
discrepancy appears to be due to two effects of the wake 
dynamics ori the measured chemistry: (1) the effect of 
averaging over species fluctuations on a scale smaller 
than can be resolved by the sampling rate, and (2) large 
turbulent mixing rates near the plume edges dramati- 
cally modifying the local equilibrium chemistry predic- 
tions. Both are accentuated by the dramatic difference 
in HO 2 concentration inside and outside of the wake 
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plume. 

We have discussed the fluctuation problem princi- 
pally from the point of view of a sampling size diffi- 
culty, but it is clear that any wake model has the same 
problem in averaging over large fluctuations within its 
representative volume elements. For a plume model 
with volume elements consistent with the data sam- 
ples, we have confirmed the speculation of Tan et al. 
[1998, p. 1723] that “turbulent fluctuations in the 
plume may significantly reduce reaction rates between 
species which are anti-correlated in distribution, such 
as between NO and HO2.” Even in our unsteady. 3-D, 
high-resolution model, there are probably some volume 
elements around the edges of the wake where fluctua- 
tions with significant anticorrelation between NO and 
IIO2 occur. To correct this, it would be necessary to use 
even higher resolution, or to introduce the major effect 
of the fluctuations into a subgrid model of the chemical 
kinetics. Since the reaction in (6) depends on the aver- 
age of the product of NO and HO2 rather than the prod- 
uct of their averages, the effective reaction rate will be 
increased (reduced) by fluctuations when they are corre- 
lated (anticorrelated). Progress has been made toward 
generating models which include the influence of the 
unresolved covariance of the reactants within a model 
[e.g., Sykes et al 1994], but a robust model for gen- 
eral conditions remains a challenging problem. Tn our 
high-resolution model we believe neglecting such effects 
only influences the sharpness of the wake edges, but in 
2-D models where ail of the downstream 3-D structure 
is averaged out in the model wake representation [e.g., 
Anderson et al 1996], the uncertainty expected should 
be similar to that found in the in situ measurements. 

In contrast to the strong influence the fluctuations 
have on HO* chemistry, the influence on the NO* - O3 
chemistry is relatively benign. For the conditions of this 
simulation, dilution of the wake tends to be faster than 
the rate at which NO within the wake destroys the 0 3 , 
so a relatively small O3 deficit is generated within the 
wake. It is not sufficiently strong to yield a strong anti- 
correlation between NO and O3 within the fluctuations. 
Thus the basic effect of the turbulence is straightfor- 
ward dilution of the wake which can be captured by 
relatively simple Gaussian plume models of the wake 
[e.g., Karol et al.. 1997]. This would not be true if the 
emission of the NO* was an order of magnitude larger 
or the dilution rate an order of magnitude smaller. 

The simulation presented was performed on a desk- 
top workstation (a 566 MHz Alpha CPU with 512 Mb 
of ram) using just over 4 days of CPU time, so it is quite 
feasible to do other simulations dealing with more com- 
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plicated chemical pathways of possible interest, and/or 
eariy jet or late plume chemistry. The transition regime 
between plume models and 3-D global models presents 
a particular challenge given the potential importance of 
turbulent mixing and large species fluctuations occur- 
ring on scales below typical grid lengths, of order 100 
km, which can be employed. We would not recommend 
attempting such a simulation at resolutions capable of 
resolving individual plumes while maintaining around 
100 possible chemical interactions, but it should be pos- 
sible to include the chemical kinetics associated with 
10-20 critical species to estimate the uncertainty asso- 
ciated with neglecting all of the subgrid fluctuations. 
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